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Abstract 


This paper presents the technical background of the Trim and Static mod- 
ule of the Matlab® Stability and Control Toolbox (MASCOT). This module 
performs a low-fidelity stability and control assessment of an aircraft model 
for a set of flight critical conditions. This is attained by determining if the 
control authority available for trim is sufficient and if the static stability char- 
acteristics are adequate. These conditions can be selected from a prescribed 
set or can be specified to meet particular requirements. The prescribed set 
of conditions includes horizontal flight, take-off rotation, landing flare, steady 
roll, steady turn and pull-up/push-over flight, for which several operating con- 
ditions can be specified. A mathematical model was developed allowing for 
six-dimensional trim, adjustable inertial properties, asymmetric vehicle lay- 
outs, arbitrary number of engines, multi- axial thrust vectoring, engine (s)-out 
conditions, crosswind and gyroscopic effects. 

Nomenclature 1 

ciij jth orientation angle of the ith engine 

B Body reference frame 

£ Earth reference frame 

T Fuselage reference frame 

h* angular momentum of the ith engine 

h ji j-terrn of the angular momentum of the ith engine 

/ Identity matrix 

/ Inertia tensor 

L ac Transformation from frame C to frame A 
n Load factor 

Nj Normal force at j'tli landing gear 
0 Nominal center of gravity 

Vi ith Propulsion reference frame 
r Residue 

v m/n Position vector of m relative to n 
s Trimmed state 

s State vector 

V Magnitude of the aircraft’s velocity relative to the wind 

V Target speed relative to ground 

W Wind reference frame 

5 Control surface deflections 

7 Target angle of climb 

yU Friction coefficient tire-runway 

1 See reference [1] for additional symbols. 
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v Wind velocity vector 

u Angular velocity vector 

uj Desired magnitude of the angular velocity (scalar) 

cr Target course angle 

r Bank angle 

£ Geometrical variables 

(i) Referring to the ith component of a vector 

Subscript 

a Referring to an aerodynamic term 

CG Center of gravity 

g j Referring to the jth contact force 

i Referring to a inertial term 

pi Referring to the ith engine 

w Weight 


1 Introduction 

Aircraft stability and control requirements are crucial for safety and perfor- 
mance. These requirements, however, are evaluated quantitatively only in the 
later stages of the vehicle’s design cycle. Unfortunately, at that stage, the basic 
aerodynamic characteristics are costly to change, resulting in an aircraft that 
is either overly conservative or that requires substantial stability and control 
augmentation to correct deficiencies. The flight control system is commonly 
expected to rectify, if possible, the legacy of stability and control deficiencies 
left by the aircraft designer. Several studies [2-4] on stability and control 
during conceptual design are available. However, most of the methods are 
ad-hoc or limited in their generality, extendability, and scope. Literature on 
high-fidelity analysis is plentiful, but the methods can not easily be adapted 
to the early phases of the vehicle’s design. 

The Matlab® Stability and Control Toolbox (MASCOT) is being devel- 
oped by the National Institute of Aerospace and the Dynamic Systems and 
Control Branch of NASA Langley Research Center. MASCOT integrates sta- 
bility and control considerations into a multi-disciplinary conceptual design 
framework. In addition, MASCOT has been conceived and developed into a 
tool that is well suited to assess stability and control characteristics at pro- 
gressive fidelity levels during the design cycle. Studies on the aircraft’s con- 
trol authority and flying qualities are made in a variable-fidelity environment, 
where the refinement of the models and of the resulting assessments progresses 
as the design cycle evolves from the conceptual to the detailed design phase. 
MASCOT is organized into five interconnected modules. The Trim and Static 
(TS) Module generates a low-fidelity assessment for a set of Flight Critical 
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Conditions. The TS Module determines if the control authority available for 
trim is sufficient and if the static stability characteristics are acceptable. The 
Open-Loop Dynamics (OD) module, generates a medium-fidelity assessment 
where a dynamic analysis of the linear and non-linear open-loop system is car- 
ried out. This includes time-based simulations and sensitivity analyses. The 
Closed-Loop Dynamics (CD) module generates a high-fidelity assessment of 
the linear and non-linear closed-loop system. This implies the integration of 
user-designed controllers, the ability to do conventional control studies and 
the assessment of the vehicle flying qualities. The Sizing and Control Alloca- 
tion (SA) module integrates MASCOT into a multi-disciplinary vehicle’s sizing 
framework. This allows the sizing of control surfaces, engines and horizontal 
and vertical stabilizers. The SA module also supports performance and ma- 
neuverability analyses and control allocation studies. Finally, the Uncertainty 
Management (UM) module accounts for the effects of low- and medium-fidelity 
input data. This paper describes only the TS module. 

This paper is organized as follows. Section 2 presents the fundamental 
mathematical model and implementation details used by all the modules of 
MASCOT. Section 3 details the trim analysis procedure and section 4 presents 
the static stability analysis capability. Finally, some concluding remarks are 
provided in section 5. 


2 Aircraft Dynamics 

Derivation of the equations of motion from first principles is available in mul- 
tiple references [1,5]. This section presents relevant information, including 
some mathematical preliminaries, on the manner in which the equations are 
implemented. Omitted definitions and developments can be easily found in 
the literature. 

2.1 Frames of Reference 

Let x- 4 , y- 4 , and z A denote the members of the orthonormal basis of frame A. 
The main frames to be used are: 

2.1.1 Flat Earth (£) 

Frame fixed to Earth where x. £ and y £ span the horizontal plane while z £ 
points downward. This frame is assumed to be inertial. 

2.1.2 Fuselage (JF) 

Frame fixed to the fuselage with origin at the arbitrary point O fixed relative to 
the fuselage, y^ aligned with the right wing and x^ aligned with the forward 
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longitudinal direction of the vehicle. 

2.1.3 Body ( B ) 

Frame fixed to the fuselage with origin at the center of gravity CG, with the 
same orientation as T . 


2.1.4 Stability (S) 

Frame fixed to the fuselage with origin at O 2 and x 5 in opposition to the 
projection, onto the x B z B plane, of the wind velocity relative to the vehicle. 
z s points downward when x s is in the horizontal plane. 

2.1.5 Wind (W) 

Frame fixed to the fuselage with origin at O and x vv in opposition to the 
velocity of the wind relative to the vehicle. z vv points downward when x vv is 
in the horizontal plane. 


2.1.6 Propulsion of the ith engine (V t ) 

Frame fixed to the ith engine having x Pi aligned with its direction of thrust. 
z Vl points downward when x Pi is in the horizontal plane. 


2.2 Transformations Between Reference Frames 


Rotation matrices of direction cosines are defined as 


LJ-) t 


1 0 0 
0 cos(-) sin(-) 

0 — sin(-) cos(-) 


Ly(') = 


cosf 

0 

sinf 


0 — sin(- 

1 0 

0 cos(-) 


LA-) = 


cosf- 
- sin( 
0 


sin(-) 0 
cos(-) 0 
0 1 


( 1 ) 

( 2 ) 

( 3 ) 


Relations among frames of reference are shown in Figure 1. The arguments 

2 The origin O is commonly chosen to be at the CG, making the frames T and B coincide. 
This is a particular case of the one considered here. The framework chosen allows performing 
studies where the location of the CG varies, without altering the procedure by which the 
aerodynamic forces and torques are calculated. 


3 




Aerodynamic 
forces and 
torques 


Forces and torques 
from the ith engine 


Figure 1. Frames of reference. Comments refer to quantities usually described 
in such a frame. 


of the rotation matrices are the Euler angles, ift, 9, 0; the angle of attack 
a , the sideslip angle /3, and the thrust vectoring angles of the ith engine 
an and a i2 . The matrix transformation between two frames can be formed 
by multiplying the transformation matrices corresponding to the reverse of 
the ordered sequence of rotations that relate the frames [5]. For instance, 
the transformation from the Wind frame W to the Earth frame S is given 
by L £w = = L z (— r ilj)Ly(—6)L x (—(j))Ly(a)L z (—/3). The reverse 

transformation is given by the inverse matrix, e.g. L Wi = (L fVV ) -1 . Note that 
since a translation relates frames T and B, no rotation matrix is indicated. 
Figure 2 shows the frames’ orientation with respect to each other. 


2.3 Preliminaries 

Constituents of the equations of motion are briefly introduced next. 


2.3.1 Aerodynamics 

Let k B be a component of a force or torque applied to the vehicle by the air. 
A linear aero-model, based on small perturbation theory, uses 


k B 


x - dk B 


( 4 ) 
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Figure 2. Reference Frames. 


where k B is a dimensional aero-coefficient, Xi is a state variable, dk B /dx% is 
a dimensional stability or control derivative, and A x* is the perturbation of 
the ith state variable from its value at trim. Evaluating this expression re- 
quires the transformation of coefficients and derivatives from >V to B and their 
dimensionalization [1,5]. Low- fidelity aero- models usually retain a subset of 
the states in Equation (4). For instance, while longitudinal forces/torques are 
usually assumed to be dependent on u, w (or a), q, w (or a), and 5 e (eleva- 
tor deflection) only, lateral dynamics are usually assumed to be dependent on 
v (or (3 ), p, r, 5 a (aileron deflection), and 5 r (rudder deflection). A general 
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aerodynamic model takes the form 


k B (h, V, a, 0, p, q, r, V, a, fi, p, q, r, 5 , 6, £, £) 


where h is the altitude, 5 contains the variables that parameterize the control 
surface deflections and propulsion settings, and £ contains the variables that 
parameterize geometrical changes in the vehicle’s layout, e.g. landing gear 
retracted or deployed. The aerodynamic forces and torques in B are given by 




_i_ r yv 

' T 0/cg 


X f ; 


w 


(5) 


where r Q/ cg is the position vector of point O relative to the CG. 


2.3.2 Weight 

The force and torque resulting from gravity are 

f J = L B£ [0,0 ,mg] T 

tw B = o (6) 

where m is the vehicle’s mass and g is the gravitational constant. 


2.3.3 Propulsion 

The force and torque generated by the ith engine are given by 

f B _ t BVif 
A p i -Lj A pz 


t p — t 
u pz L pz 


+ V 0/cg X fpi 


(7) 


where 


t t _ j TV, 


t Ti y'Pi X f . 
L P* ' 1 e;/0 A 1 P* 


Vi 


and r^ l / 0 is the position vector of the ith engine relative to point O. 

2.3.4 Contact forces 

The force and torque generated by the ground on the jth landing gear are 
given by 


f B = f T 


f B _ (JB , 

— 1 r n/m + r 


O/cg 1“ ‘-cpj/O^j 


X 7 g j 


( 8 ) 


where r ^ Pj / 0 is the position vector of the contact point between the ground 
and the jth landing gear relative to point O. 
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2.3.5 Kinematics 


The relation between the Enler angles and the angular velocity is given by 


w® = 

’ p " 

q 

— 

' 1 

0 

0 

COS (0) 

— sin(d) 
sin (0) cos (9) 


" 0 " 

6 


r 


0 

— sin(0) 

cos(^) cos(0) 


. . 


2.3.6 Inertia 

Denote by w B g = [u, v,w] T + L B£ v £ the vehicle’s velocity relative to ground, 
where z/ , the wind velocity, is assumed to be constant. The inertial force is 
given by 

f ; B = m (a® + u B x v® ) (10) 

where a.® = [u, v, w] T + L B£ v £ is the acceleration of the vehicle relative to 
the ground in frame B, and uj b is the angular velocity of frame B relative to 
£. Notice that Coriolis and the relative acceleration terms are neglected due 
to the flat-earth assumption. The velocity and acceleration can be written in 
terms of the wind-relative speed V, the angle of attack a, the side-slip angle 
/3, and their time derivatives according to 

[u, v, w} 1 = V [cos(a) cos (/3), sin(/3), sin(a) cos(/3)] T 

The inertial torque is given by 



where 1 B is the inertia tensor, u B = [p, q,r ] 7 is the angular velocity, Co B = 
[p, q,r] T is the angular acceleration and h,; K is the angular momentum of the 
ith engine. These terms can be calculated from 

= L FPi Y\ 1 j Pi + h 2i JP 

I s = F + m [(r B cg/0 ■ r B cg/0 )I - r B g/ 0 r B g/0 T ] 

where I is the identity matrix while hi j Pl and h 2 ^ are terms used to describe 
the angular momentum of the rotating components of the ith engine. Usually, 
one of these two vectors is zero. If the axis of rotation of the engine moves 
relative to the fuselage h 2 i~ F = 0, e.g. rotating nacelles; otherwise = 0, 
e.g. fixed nacelles. 


2.3.7 Navigation 

The navigation equation is given by 


( 12 ) 
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2.4 Equations of Motion 

Define the target velocity vector as 

v 5 = L*(-(7)L y (-7) [F, 0,0] T 

where V is the target speed relative to ground, a is the desired course angle 
and 7 is the desired angle of climb. Figure 2 shows the relations among frames 
and relevant quantities. The force, torque and navigation equations can be 
written as 


D = fa S + U B + Y J f P^ + ( 13 ) 

i 3 



(15) 


where specific references to vector components are given in Equation (15), 
e.g. vf (1) is the first component of vf . We will refer to for k — 1, 2, 3; as 
residues. Providing that Equation (9) holds, the force and torque equations are 
satisfied when iq = 0 and r 2 = 0 respectively. The vehicle’s velocity will match 
the target velocity when r 3 = 0. The complete set of differential equations 
that describe the vehicle’s dynamics is given by Equations (9,12,13,14). 


3 Trim Analysis 

Trim analysis requires searching for the states that satisfy a simplified realiza- 
tion of the equations of motion. This search is to be done for all flight critical 
conditions. A flight critical condition is defined as the combination of a flight 
condition and an operating condition. Flight sets of interest are introduced 
next. 

3.1 Flight Conditions 

A flight condition is an instantaneous flight maneuver prescribed by Equations 
(13,14,15) with r*. = 0 for k = 1,2,3; along with some kinematic and kinetic 
constraints. Physical interpretations of flight conditions of interest and the 
corresponding constraints are presented next. 



3.1.1 Horizontal Flight 

This flight condition occurs when the vehicle’s velocity relative to ground is 
horizontal. This implies y = 0 while lo b , Cj b , fg^ for all j, and a B g are zero. 

3.1.2 Take-Off Rotation 

This flight condition occurs when the vehicle, moving at lift-off velocity on the 
ground, starts rotating about the axis of rotation of the main landing gear. 
At that instant, the normal forces at the wheels whose axes do not coincide 
with the axis of rotation are zero. This condition implies 6 = 0, 7 = 0, 0 > 0 
and c o B = 0. In addition, the constraints fj < /oNj, where fj is the friction 
force, p is the friction coefficient between the tire and the runway and Nj is 
the normal force; bound the components of fg^. 

3.1.3 Landing Flare 

This flight condition occurs when the vehicle, flying at approach speed, starts 
a rotation that raises its nose. This is done to prevent touching down at high 
sink rates. This condition implies 6 > 0, uo B = 0 and fgj^ = 0 for all j landing 
gear. 

3.1.4 Steady Roll 

This flight condition occurs when the vehicle rolls at constant angular velocity. 
This implies u B = [p, 0,0] T while u B , f gJ A for all j and a B g are zero. Another 
roll condition of interest is given by u B = L BW [Co, 0, 0] T , case in which the axis 
of rotation coincides with the direction of the aircraft’s velocity relative to the 
wind. 

3.1.5 Steady Turn 

This flight condition occurs when the vehicle turns at a constant, vertical 
angular velocity. This implies u B = L B£ [0,0,u;] T , v^(3) = 0, 7 = 0, a B g = 
L B£ (lo £ x vf 9 ), Co B = 0 and = 0 for all j. If the resultant of the weight 
and the inertial force is in the x H z B plane, condition commonly refereed to 
as coordinated turn, the approximation u B = ta a{r)g/V \—6,n cos(0), l/n] T , 
where r is the bank angle [6], and n is the load factor; is commonly assumed 
for v £ = 0. Note that wind velocity does not contribute to the centripetal 
acceleration of the turn. 

3.1.6 Pull-Up/Push-Over 

This flight condition occurs when the vehicle flies at the bottom or the top 
of a curved path and the centripetal acceleration remains vertical. These 
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conditions imply 7 = 0, vf (3) = 0 and 

lu b = (n [ 0 ,- 1 , 0 ] T 

while fgj^ = 0 for all j and a B g = L B£ [0, 0, (n — 1 )g] T ■ Non-maneuvering, i.e. 
n = 1, pull-up, i.e. n > 1, push-over, i.e. n < 1, and ballistic, i.e. n = 0, 
flight conditions are special cases of interest. Note that wind velocity does not 
contribute to the centripetal acceleration. 

3.2 Operating Conditions 

An operating condition determines the circumstances occurring during a given 
flight condition. Engine(s) out, crosswind, and crabbed configurations exem- 
plify this aspect. This is attained by means of additional constraints, e.g. an 
engine out condition requires = 0, f t f = 0, h-xf 1 = 0 and h = 0 for 
a given i\ and parameter assignments, e.g. crosswind requires v £ ^ 0. We will 
refer to quasi-static FCC when the intended acceleration is non-zero. 

3.3 State Partition and Trim Search 

Denote by s a vector containing all the state variables. The state variables 
are the Euler angles {0,0,0}, the velocity components {V, a,/?} 3 , the control 
surface deflections {5}, the angular velocity components {p, q,r}, the propul- 
sion setting {an, d 2 i,f P i ' F for all engines, the contact force {fg^} at all 
landing gears, and the quasi-static terms 

Equations (13,14,15) define a coupled system of nine scalar equations. For 
a given FCC, the vehicle will be trimmed when there exist a realization of 
the state, denoted by s, that leads to r*,(s) = 0 for k = 1,2,3. The trim 
conditions found using the TS module can be used as reference conditions for 
the sizing and control allocation module (SA) for both low and progressive 
fidelity models. This will allow for the evaluation of critical maneuverability 
requirements. 

The search for the trimmed states 's is done by minimizing the sum of the 
norms of the residues. When the number of states exceeds nine variables, the 
resulting system of equations is under-determined. To prevent solving such a 
system and to pose flight critical conditions that impose stringent demands 
on the control solution, a partition of the state vector is required. Hence, the 
state vector is partitioned into three groups: the group of unknowns, the group 
of constants and the group of related variables. The latter group contains 
all the state variables assuming a value equal to the value of an unknown. 
There is an equality associated to each member of this group. Each flight 

3 Equivalently, the aircraft’s velocity relative to the wind can be prescribed by {u, v, re} 
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condition set, with the exception of the user defined case, has a prescribed 
state partition. Special care should be given when s is partitioned by the User 
since the resulting FCC might be physically meaningless. 


4 Static Stability Analysis 

The Trim and Static (TS) module performs a static stability analysis for each 
FCC, by calculating equivalent stiffnesses, dampings, and static margins. The 
signs of the stiffnesses determine the direction of the initial restoring force 
from a perturbed state. The sign of the damping terms indicate if energy 
is flowing either from or to the vehicle. The static margin is a metric that 
assess the CG location in pitch by measuring the closeness to neutral stability. 
Notice that a static stability analysis refers to the instantaneous response of 
the vehicle from a perturbed state in the vicinity of s and does not constitute 
a stability analysis in the dynamic sense. Statically unstable vehicles could 
be dynamically stable and vice versa. Details on the static metrics are given 
below. 


4.1 Equivalent Stiffnesses 

Static stability in roll, pitch and yaw are attained when the following inequal- 
ities are satisfied 


/ *9r 2 (l) 

V <90 

/ <9r 2 (2) 
\ da 
( <9r 2 (3) 

V dp 


< o 

< 0 
< 0 


4.2 Equivalent Dampings 

Static stability in roll, pitch and yaw requires 


dCl 


dp J 8=1 

/ dCm\ 

vw)^ 

( a -p) 

V ,Jr 


< 0 


< 0 


< 0 


(16) 

( 17 ) 

(18) 


(19) 

( 20 ) 
( 21 ) 
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where Cl, Cm, and Cn are the non-dimensional aerodynamic coefficients for 
rolling, pitching and yawing respectively, and p, q and r are the non dimen- 
sional components of the angular velocity. Notice that these dynamic deriva- 
tives are required to build the linear aerodynamic model in Equation (4). 

4.3 Longitudinal Static Margin 

The static margin is determined by the longitudinal distance between the CG 
and the vehicle’s neutral point. The neutral point is defined as the location 
where the equivalent pitch stiffness is zero, i.e. where Equation (17) no longer 
holds. The static margin is provided as a percentage of the mean aerodynamic 
chord. Usually, subsonic transport vehicles have positive margins while fighters 
have negative margins. 

5 Concluding Remarks 

This paper presents the mathematical background of the Trim and Static Mod- 
ule of MASCOT. Even though this module generates low-fidelity stability and 
control assessments, the developments presented also support the medium- 
and high-fidelity tasks performed by other modules. MASCOT assessments 
should not only evolve in fidelity as the design cycle does but more impor- 
tantly, they should be performed according to the resources and information 
available. For instance, if the conceptual designer is not willing to spend more 
than few minutes in the evaluation of a concept, longer turn-around times 
will certainly render this tool impractical. These requirements as well as the 
underlying consistency among the variable fidelity assessments constitute the 
main challenges of this endeavor. 
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are adequate. These conditions can be selected from a prescribed set or can be specified to meet particular requirements. The 
prescribed set of conditions includes horizonta flightl, take-off rotation, landing flare, steady roll, steady turn and pull-up/ 
push-over flight, for which several operating conditions can be specified. A mathematical model was developed allowing for 
six-dimensional trim, adjustable inertial properties, asymmetric vehicle layouts, arbitrary number of engines, multi-axial thrust 
vectoring, engine(s)-out conditions, crosswind and gyroscopic effects. 
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